Role of electronic structure in photoassisted transport through atomic-sized contacts 



O 

o 

(N 



c3 



i 

C 

o 
o 



> 

in 
o 
in 

o 

O 

-i— > 
o3 



i 

c 

o 
o 



J. K. Viljas 1 ' 2 and J. C. Cuevas 3 ' 1 ' 2 

1 Institut fur Theoretische Festkorperphysik, Universitat Karlsruhe, D-76128 Karlsruhe, Germany 
2 Forschungszentrum Karlsruhe, Institut fur Nanotechnologie, D-76021 Karlsruhe, Germany 
3 Departamento de Fisica Teorica de la Materia Condensada, 
Universidad Autonoma de Madrid, E-28049 Madrid, Spain 
(Dated: February 6, 2008) 

We study theoretically quantum transport through laser-irradiated metallic atomic-sized contacts. 
The radiation field is treated classically, assuming its effect to be the generation of an ac voltage 
over the contact. We derive an expression for the dc current and compute the linear conductance 
in one-atom thick contacts as a function of the ac frequency, concentrating on the role played by 
electronic structure. In particular, we present results for three materials (Al, Pt, and Au) with very 
different electronic structures. It is shown that, depending on the frequency and the metal, the 
radiation can either enhance or diminish the conductance. This can be intuitively understood in 
terms of the energy dependence of the transmission of the contacts in the absence of radiation. 

PACS numbers: 73.63.-b, 73.50.Pz, 73.63.Rt, 73.40.Jn 



I. INTRODUCTION 

The study of electronic transport in microscopic and 
nanoscale electrical contacts under the influence of time- 
dependent external fields has a long history. Perhaps as 
the most famous example, superconducting tunnel junc- 
tions subjected to microwave radiation exhibit a step-like 
structure in their current-voltage (/-V) characteristics 
This can be understood in terms of inelastic ("pho- 
toassisted") transport of electrons across the junction. In 
the early theoretical work of Tien and Gordon (TG) 0, 
this phenomenon was described by a harmonic voltage 
at the radiation frequency lu applied to one of the leads, 
giving rise to photo- "sidebands" associated with the "ab- 
sorption" or "emission" of an integer multiple of the pho- 
ton energy fko. This theory has also been extended to 
describe superconducting atomic point contacts [3[ and 
the predictions have been confirmed experimentally us- 
ing microfabricated Al break-junctions under microwave 
irradiation [J]. In addition, a similar approach has been 
used to describe laser-irradiated junctions in scanning 
tunneling microscopes (STM) H. In these, as a result of 
inherent asymmetries in the geometry and the materials 
of the junction, laser irradiation can cause dc (rectifica- 
currents even in the absence of a dc bias voltage 



Over the years, several types of model calculations 
have been employed also to describe the ac-response of 
semiconductor heterostructures [1, Q and other meso- 
scopic systems flOl. ITU Itj , as well as atomic and molecu- 
lar contacts [l3l Il4 lla. Ho. fl"^ . Among these, the TG-like 
theories have been quite successful in gaining a qualita- 
tive understanding of light-induced currents d, [l3| . 

Metallic atomic-sized wires fabricated with STM or 
break-junction techniques have turned out to be ideal 
systems for investigating electronic transport at the 
nanoscale 18]. The bulk of the research in this field 
so far has concentrated on stationary transport proper- 
ties, while systems being driven by time-varying external 



fields (such as laser light) have received less attention. 
Very recent experiments on laser-irradiated gold contacts 
support the idea that photoassisted processes may play 
an important role in their transport properties [19[ . To 
take the first steps towards a microscopic description of 
experiments of this type, we address in this paper the role 
of electronic structure in photoassisted transport through 
atomic-sized junctions. This problem is not only relevant 
for the field of atomic contacts, but also for molecular 
electronics, where the role of the metallic contacts on 
photoassisted transport through molecular junctions re- 
mains to be understood. 

In different types of metals, the nature and number of 
conduction channels in one-atom contacts reflect the va- 
lence of the metal, i.e., what type of atomic orbitals are 
available at the Fermi energy [20]. How is this difference 
between the metals seen in their response to irradiation? 
To investigate this question, we use a tight-binding model 
and explore mainly one-atom thick contact geometries 
like the one shown in Fig. [TJa). We are interested in the 
linear conductance Gdc = dI/dV\v=o as a function of the 
radiation frequency uj, which we call photoconductance. 
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FIG. 1: (Color online) (a) Two infinite fee [001] surfaces con- 
tacted with pyramids of 10 atoms in each, forming a "dimer" 
contact. The tip-to-tip distance is denoted with D, all other 
interatomic distances are the same as in bulk, (b) Two model 
ac voltage profiles £/f c , where i labels atomic sites and or- 
bitals in the C region: a double step (A) and a linear ramp 
(B) between the lead values C/£ c and Ur c . 
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Following the TG ideology, we model the effect of radia- 
tion with a time-periodic voltage, but we also take into 
account the effect of the voltage profile across the contact 
[see Fig. [IJb)]. While many previous model calculations 
with periodic driving fields are based on Floquet theory 
[H> US], our method is based on non-equilibrium Green 
functions @, @, HH H2| ■ Our approach allows for a real- 
istic description of the photoconductance of atomic con- 
tacts, with a single free parameter describing the local 
intensity of the radiation. 

As examples of metals with very different electronic 
structures we choose Al, Pt, and Au. We study one-atom 
thick junctions of these materials mainly in the so-called 
contact regime, but also in the tunneling regime. In the 
case of the contact regime we find the following results. 
For Pt the effect of the ac voltage is almost always to 
decrease Gd c from its value in the absence of radiation. 
This is because the Fermi energy ep lies at the edge of 
the d band, and exciting electrons above ep necessarily 
decreases their transmission probability, since less trans- 
mission channels are available there. For Al, where cf 
is at the beginning of the p band, also an enhancement 
of the conductance is possible. The magnitude of these 
changes depends on the intensity of the radiation, but 
can be up to 10% or more at visible frequencies. For 
Au, ep is in the s band, with a single completely open 
transmission channel for a wide range of energies. Thus, 
low-frequency radiation has no effect on the conductance, 
while in the visible range both an increase and a decrease 
are possible. These conclusions are based on detailed nu- 
merical simulations, but they can be understood in an 
appealing way in terms of the energy dependence of the 
transmission in the absence of the ac drive. 

The rest of the paper is organized as follows. In SecUTl 
we introduce the theoretical model used to describe the 
electronic structure of atomic-sized contacts and the ef- 
fect of irradiation. Section HTTl is devoted to the derivation 
of a general formula for the dc current through atomic 
contacts subjected to an ac voltage. In Sec. IIVI we de- 
scribe examples of numerical results of the frequency de- 
pendence of the linear conductance for one-atom thick 
contacts of Al, Pt, and Au. Finally, in Sec. [V] we discuss 
some of the assumptions and restrictions of our model 
and present the main conclusions. 

II. MODEL 

We describe the laser-irradiated point contact with an 
spo?-tight-binding Hamiltonian of the form 

H(t)=H + H 1 (t) 
H = J24 H ijdj, H 1 (t)=Y J d\W l3 (t)d ] , C 1 ) 

ij ij 

where the indices i,j run over the different atoms and 
orbitals, including the two spin directions. The creation 
and annihilation operators (dj and dj) satisfy {di,dj} = 



and {di,d^} — [S^ 1 ]^, where S is the overlap matrix 
of the non-orthogonal basis [23| ■ The Hamiltonian Hq is 
for the system without radiation and dc voltage, while 
Hi(t) includes them both. The matrix elements of the 
Hamiltonian and the overlap matrix ( Hq and Sij) are 
taken from the parametrization of Ref. [24]. All matrices 
in the spin-orbital basis (such as H and S) are denoted 
with a boldface symbol. 

We consider ideal symmetric geometries of the type 
shown in Fig. [lja), with a single-atom thick constriction. 
This type of "dimer" structure is suggested by molecu- 
lar dynamics simulations as the most common one in the 
last conductance plateau (25[. When the distance D be- 
tween the tip atoms corresponds to the bulk interatomic 
distance, the junction is said to be in the contact regime. 
We also study larger D values, where the junction enters 
the tunneling regime. For the calculation of transport 
we shall, in the usual way, divide the system into left 
lead (L), central (C), and right lead (R) regions. The 
leads are modeled with infinite surfaces, where the fee 
[001] axis coincides with the transport direction. We in- 
dex the regions with the label = L, C, R. As a matrix 
index, this label also indicates collectively all the orbital 
indices of the respective region [23| . 

We assume that the L and R lead potentials are spa- 
tially constant and harmonic with angular frequency to, 
such that T u = 2ir/u> is the oscillation period. Thus 
Wnn{t) = U n (t)Snn, where U a (t) = U*f + t^ c cos(wt), 
and Q = L,R, The applied dc voltage V = {U^-U^/e, 
where — e is the electron charge, is assumed to be in- 
finitesimal. For the lead-center hoppings we also as- 
sume Wcsi(t) = Un(t)Scn, etc. The central poten- 
tial is assumed to be of the form Wcc(t) = Wq C + 
W% c c cos(u;t), or [W C c(t)}ij = [S C c]ij(Ui(t) + Uj(t))/2 
[IJ. Here, U t {t) = Uf c + Uf c cos(wi) is the same for all 
orbitals i on the same atom in region C. The actual shape 
of Ui(t) within C should in principle be obtained self- 
consistently through the solution of a Poisson equation, 
such that screening and local field-enhancement effects 
would be properly accounted for. In this case the appear- 
ance of higher harmonics of uj and phase shifts in W cc {t) 
would be possible. As this would be computationally very 
costly, we shall do the following: the dc part Uf c is fixed 
by a requirement of local charge neutrality in equilibrium 
[2^,[2^| and incorporated into Hqc, while for f/" c we just 
assume some simple forms. Below we consider two model 
profiles [see Fig. [Tib)]: Uf c = (C/£ c + U%)/2 for all or- 
bitals i in region C (profile A), and one which linearly 
interpolates between Uff and U^ 2 (profile B). For the 
symmetric junctions that we are considering, symmetric 
profiles of this type are the most reasonable. 



III. CURRENT FORMULA 

We only consider the time-averaged current, and can 
thus neglect displacement contributions without losing 
current conservation [ToL [ill [l2| . The particle current is 
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given in terms of non-equilibrium Green functions by @ 
rT " dt 



I = - 



— Trl(G< L ot LC ~tcLoG< c )(t,t)}. (2) 



Here, in the case of the non-orthogonal basis [23| 

t C L(t, t') = [{H CL + W CL {t)) - ScLihd t ]h8{t - t'), 

and using the "Langreth rules" , Gq L = G r cc o t CL o 
g^L + Gp C otcL° @- The product o is defined by 
(AoB)(t,f) = J(ds/h)A(t,s)B(s,t'), where A and B 
are matrices in the spin-orbital basis, over which the trace 
Tr acts. The Green functions Gqq, with x = r, a, ^ and 
f2, Q' = L,C,R are defined as usual [9|, l23j, and <7qq are 
the functions for uncoupled leads. The CC component 
of the full retarded function satisfies 

[Seeing - Hoc - W C c(t)]G r cc (t,t') = h5(t - t') 



whereas Gp C 



-CE r L oG r cc )(t,t') + (^ R oGhc)(t,t'), 



G 



cc 



(3) 

Here we 



defined the "lead self-energies" £ = *csi ° 9no tnc, 
with fl = L,R. The solutions of 

[Snnift^ - ifnn - W Qn (t)]g r nn(t, = M(i - f ) (4) 



lead to self-energies of the form 



2tt 



-fc(t-f)/» E £««( e ). (5) 



The equilibrium self-energies are given by S^ e, (e) = 
<cn(e)9an 9 ( e )*^C'(e), where the Green functions g^ni 6 ) 
of the infinite surfaces are obtained by a decimation 
method [23[, while icn(e) = i?cn - eScn and tnc = 
[tcnV- We also have g^(e) = -/(e) [ff^(e)-ff^(c)] 
where /(e) is the Fermi function. 

Since the Green functions and self-energies A(t,t') are 
all periodic in the time (t + t')/2, we can simplify the 
analysis by working in Fourier coordinates. We define 
the harmonic matrices A(e) with components L4] m , n (e) = 
A m _„(e + (to + n)fkt>/2), where to, n are integers, and 



A n (e) 



dT 



inuT 



dr 



A(T + r/2, T - r/2). 



J h 

In this way we find = G r (X]f + £%)G a , while 



[G r 



\eS, 



cc 



H cc ]l-W-± 



r,a 
L 



(6) 



Here [e] m ,„ = (e + mtkj)5 m>n , [W] m , n = W a f c (S rn - 1 ^- l + 
3m+i.n)/2, and we have dropped the indices CC from 
the harmonic matrices. The self-energies as well as the 
scattering rates To = i(^Q — Y,q) are related to the cor- 
responding equilibrium quantities by 



[T n 



where we define the components 
[f^L,n(e) = Jm+z (an) J n +i (an) r£(e - - Hiw), 



with a similar equation for E^^(e). Here Tj 
i[E r ' eq (e) - £ Q ' e<z (e)], J; are Bessel functions of the first 
kind, and an = U^f/hu. All of the harmonic matrices 
satisfy the symmetry [A^] m+k . n+k (e) = [A i - l+k ^] m . n (e + 
khw). Finally, we may note that J(dt/T u! )Tr[(A o 

B){t,t)] = J ^(de/2^)Tr w [i(e)B(e)], where Tr u [A(e)] = 
J2 m Tr[A]„ liTO (e). For numerical calculations, the matri- 
ces A(e) must be truncated to a few lowest indices m, n, 
but the results converge rapidly with the cutoff. 

Using the above definitions and general symmetries like 
G r ~G a =G> -G<, Eq. © yields the dc current 



^) 



I = 



/ ^E^t^G^K/f -/ 



R )■> 



(8) 



where /^(e) = /(e — U{[ c — khuj). In the absence of 
an ac field this reduces to the standard Landauer-type 
formula [27|]. Although we are not computing Ui(t) self- 
consistently, Eq. ((Sj) is still gauge invariant in the sense 
that a spatially constant potential added everywhere has 
no effect. Thus, the results only depend on Ul (t) — Un(t) 

El. 



IV. RESULTS 

The experimentally accessible quantity which we cal- 
culate is the linear conductance Gd c (u) = dI/dV\v=o- 
In what follows we shall assume zero temperature. If the 
lead self-energies Yf£ e R are furthermore assumed to be 
energy-independent ( "wide-band" approximation) , then 
the full result for potential profile A simplifies to 11 1 



G dc (uj) = G J2{Ma/2)] 2 T eq (e F + Ihiu), 



(9) 



where Go = 2e 2 /h and T eq (e) is the equilibrium trans- 
mission function [28[ . We have also defined the parameter 
a = ccl — ccr = (Ul c — Uft C )/hoj, which measures the local 
intensity of the radiation 0] . Equation §§§ describes elec- 
trons incident from the different sidebands (after having 
"absorbed" or "emitted" I photons) being transmitted 
elastically through the constriction, which mostly deter- 
mines T eq (e). In reality, the electric field is nonzero only 
in the constriction and thus the actual physical transi- 
tions must occur there. Note that Eq. ([9]) should repro- 
duce the results for profile A only in the limit of small uj 
and a. Still, it works surprisingly well for all the cases 
presented below. 

Let us first discuss the results for the contact regime, 
which is the subject of our main interest in this pa- 
per. In Fig. [5] we show results for a dimer Al contact, 
where ep lies in the 3p band. The upper panel shows 
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FIG. 2: (Color online) Top panel: Equilibrium transmission 
T eq and its decomposition into conduction channels 21,2,3,4 
for an Al dimer contact. The position of ef is indicated by a 
vertical dotted line. Lower panels: Zero-temperature photo- 
conductance for several values of a as a function of frequency 
ui using the voltage profile A (a), profile B (b) and Eq. Q 
(c). In (b) the wavelengths A with a tick spacing of 400 nm 
are shown. The range of visible light is indicated by vertical 
dotted lines. 



the transmission function T eq (t) and the lower panels 
show Gd c (w) for the two voltage profiles as well Eq. © 
for several values of a. In the absence of an ac voltage 
(u> = 0) the conductance is close to Go, and is domi- 
nated by three conduction channels due to the contribu- 
tion of 3s and 3p orbitals [2(J. At finite u>, the relative 
change 5G dc {uj) = [G dc {uj) - G dc (w = 0)]/G dc (u = 0) 
is initially negative, but can then rise to positive values 
of w 10% towards visible frequencies. This behavior is 
similar for both profiles A and B, as well as for Eq. ^j. 
From the latter result, we can interpret our findings in the 
following appealing way. For a < 1 only the first side- 
bands (I — 0, ±1) contribute to the transport. In this 
limit, according to Eq. ((9|), 8Gdc{u) measures, roughly 
speaking, the "second derivative" of T eq (e) on the scale 
of Hlu around tp. Thus, for instance, the conductance 
enhancement in the visible range follows from the trans- 
mission increase for electrons promoted above t f (due to 
"absorption" ) overcoming the corresponding decrease for 
electrons moved below tp (due to "emission"). 

Figure[3]shows the corresponding results for Pt. In the 
absence of radiation the conductance is close to 2.1Go due 
to the contributions of mainly three conduction channels, 
which originate from the 6s and 5d orbitals. In this case, 
and in the contact regime in general for Pt, the effect 
of the radiation is almost always a significant reduction 
in conductance. This is understandable, since tp lies at 
the edge of the d band, and photon absorption leads to 
an energy region where fewer open transmission channels 
are available and T eq is smaller. Note that for low u>, the 
full results are again well described by Eq. ([9]). Let us 
remark that we only compute Gdciu) for low enough lj 
and a, so that the electric fields of the radiation remain 
reasonable (< 3 ■ 10 9 V/m). 



FIG. 3: (Color online) Same as Fig. H but for Pt. 
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FIG. 4: (Color online) Same as Fig. [2] but for Au. 



The results for Au are shown in Fig. [4] As can be seen 
in the upper panel, the conductance for lo = is equal to 
IGo with a single open channel arising from the contri- 
bution of the 6s orbitals [20j. Moreover, notice that the 
transmission around tp is very flat. Due to this flatness, 
for frequencies up to fuv w 1 eV the effect of radiation is 
practically negligible. In the red part of the visible range 
(hui < 2 eV) we find that 5Gd c (u) > up to a few per- 
cent, although this depends rather strongly on the choice 
of the voltage profile. This increase in the conductance is 
due to a contribution of the 5d bands located 2 eV below 
ep, where the number of open transmission channels is 
higher than at tp. At higher frequencies SGdc{uj) < 0, 
as for Pt. We have also studied Au contacts with atomic 
chains of varying length and the results remain qualita- 
tively similar, although in the case of profile B the am- 
plitude 5Gdc(w) < becomes smaller as the number of 
chain atoms increases. 

When the distance D between the tip atoms of the 
dimer contact (Fig. [1]) is increased, the tunneling regime 
is approached. Here Gdc{^) decreases exponentially with 
increasing D, but SGd c {^) has a tendency to saturate. 
This is easy to understand, because it may be shown 
that for very large D the magnitude of the conductance 
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FIG. 5: (Color online) Tunneling limit for Al contacts. Top 
left panel: zero- frequency conductance Gd c (u) — 0) as a func- 
tion of tip distance D. The vertical lines indicate the distances 
where the examples with corresponding linestyles in the other 
panels are computed. Top right panel: transmission T eq (e) 
for the example distances. Lower panels: the left-hand (a,d), 
central (b,e), and right-hand (c,d) panels are for profile A, 
profile B, and Eq. ©, respectively. The panels (a)-(c) are for 
a = 0.5 and (d)-(f) for a = 1.0. 



is approximately determined by the square of the slowest- 
decaying hopping integral between the tip atoms, which 
enters the conductance formula as a prefactor. The form 
of <5Gdc(w) can, however, be very different from the con- 
tact regime. 

An example of the tunneling regime results for Al is 
presented in Fig. The top panels illustrate the expo- 
nential decay of Gd c (w = 0) and T eq {e) with D, while 
in the lower panels the quantity Gdc{u)/Gdc(u = 0) is 
shown for two values of a and for several distances D. 
The results again look quite similar for the two voltage 
profiles as well as for Eq. © . The quantity SGdc can ob- 
tain both positive and negative values, and its saturation 
with D is clearly visible. For Pt (see Fig. [6]) we find that 
Gd c (w) is otherwise flat, but there is a sharp resonance 
at tiw » 1 eV where 5Gd c (w) can take on positive values 
of up to a few hundred percent. This is due to resonant 
transmission through a level formed by the d orbitals 
of the tip atoms, which can be seen also in the T eq (e) 
curves. In the case of Au there exists a rather similar, 
but broader positive peak covering the visible range (see 
Fig. [7|). For each metal, we only consider large enough 
D to see the saturation of 6G^ c (u>). Indeed, the tight- 
binding parametrization we are using is based on bulk 
calculations [HI and is probably not good for very large 
interatomic distances. Furthermore, the charge neutral- 
ity shifts mentioned in Sec. [II] are strongest for the tip 
atoms. Therefore the peaks in Gd, c (ui) observed for Pt 
and Au, for example, should be taken with some reserva- 
tions. For very large D, we would also not expect to see 
such a good agreement between the results for profiles A 



FIG. 6: (Color online) Same as Fig. 02 but for Pt. 
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FIG. 7: (Color online) Same as Fig. O but for Au. 



and B, because the tunneling conductance for profile B 
depends on the local densities of states of the tip atoms 
in a way that cannot be written in the form ©. Nev- 
ertheless, our results serve as good illustrations of the 
different phenomena that may potentially arise in the 
tunneling regime. 

Above we have only shown some example results for 
very idealized symmetric geometries. In general, the 
signs and magnitudes of 8Gdc{u) depend sensitively on 
details of the atomic structure, as does T eq (e). For a 
more complete analysis one should therefore study also 
larger contacts and carry out a statistical exploration of 
geometrical variations along the lines of Ref. [2f| Based 
on Eq. © , we can still expect that in the limit of several- 
atoms-wide contacts the relative effect of the ac voltage 
gradually becomes smaller as T eq (e) becomes smoother. 
A direct comparison with the ongoing experiments [l9j 
will be postponed for later. 
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V. DISCUSSION AND CONCLUSIONS 

Some of the assumptions of our model and the effects 
not taken into account are worth discussing. First, we 
assume a flat potential in the leads, which requires a 
complete screening of the electric field. This is well sat- 
isfied in metals for frequencies u) much below the plasma 
frequency but, as uj begins to approach the visible range, 
the screening is weakened. On the other hand, as we 
have mentioned, the screening in the central region is not 
treated self-consistently. One of the main concerns here 
is that local surface-plasmon modes in small geometries 
tend to have their frequencies close to the visible range, 
and their excitation can lead to huge field-enhancement 
effects [5] . Although this is not a problem for our model, 
one should bear in mind that a is not simply related to 
the field intensity away from the contact and, following 
the spirit of the experiments in superconducting contacts 
[H, it must be understood as an adjustable parameter. 
Also, at visible frequencies "multiphoton" processes can 
already cause a photoelectric effect, i.e., excite electrons 
above the vacuum level, which typically lies 4-6 eV above 
e F . 

Heating accompanies any possible effect arising from 
absorption of light, which in metals is specially pro- 
nounced in the optical range due to the onset of inter- 
band transitions. While we may expect that the effect 
of temperature is just to broaden our results, in practice 
thermal expansion can play an important role. This is 
well documented in the STM context ||, where expan- 
sion typically brings the tip closer to the sample, thus re- 
ducing the tunneling gap width. This results in a strong 
enhancement of the tunneling current. In the case of 
atomic wires (i.e., in the contact regime), it is not obvi- 
ous in which sense and to what degree thermal expansion 
affects the conductance. Assuming that the expansion 



simply mimics a mechanical closing process in the STM 
or break-junction experiments, it can lead to either an 
increase or a decrease of the conductance, depending on 
the metal. For instance, for Al contacts, which exhibit 
raising plateaus upon stretching [291 ], one would expect 
a decrease of the conductance due to thermal expansion, 
as opposed to the effect of the electronic structure in the 
visible range (see Fig. [5]). In this sense, our predictions 
can be valuable for distinguishing in an experiment (such 
as Ref. ITii ) between the contributions of the different ef- 
fects to the photoconductance. 

In conclusion, we have modeled electronic transport in 
atomic point contacts subjected to external electromag- 
netic radiation. The radiation has been described by an 
ac voltage over the contact. Within a non-equilibrium 
Green function method, we have derived a formula for 
the dc current in the presence of such an ac drive. Us- 
ing a tight-binding model, we have applied the method 
for describing atomic-sized contacts of Al, Pt, and Au, 
and have found that the qualitative modification of the 
dc conductance by the ac voltage can be predicted from 
the equilibrium transmission function. Depending on the 
metal, the detailed structure of the contact, and the ex- 
ternal frequency, the effect can be either an increase or 
a decrease in the conductance. At present, experiments 
are under way to test these predictions [l9j . 
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